Conditional preparation of arbitrary atomic Dicke states 
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We propose an experimentally accessible procedure for conditional preparation of highly non- 
classical states of collective spin of an atomic ensemble. The quantum state engineering is based on 
a combination of QND interaction between atoms and light previously prepared in a non-Gaussian 
state using photon subtraction from squeezed vacuum beam, homodyne detection on the output light 
beam, and a coherent displacement of atomic state. The procedure is capable of non-deterministic 
preparation of a wide class of superpositions of atomic Dicke states. We present several techniques 
to optimize the performance of the protocol and maximize the trade-off between fidelity of prepared 
state and success probability of the scheme. 
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I. INTRODUCTION 

Laws of quantum physics enable to process and trans- 
mit information in ways that would be impossible or 
very difficult to achieve classically. Prime example is 
the unconditionally secure quantum key distribution that 
is already approaching the stage of commercial applica- 
tions [1]. Most of the developed quantum communica- 
tion protocols employ light as a carrier and processing 
medium. Quantum information needs however not only 
to be transmitted and processed, but it requires also stor- 
age facilities. Storage of quantum information can not 
be accomplished using classical means, but it requires a 
special - quantum - memory. Quantum memory is in- 
dispensable for construction of quantum repeaters [E El 
that combine it with entanglement distillation and swap- 
ping to efficiently distribute entanglement over lossy and 
noisy quantum channels. A very promising medium for 
quantum memory are ensembles of atoms trapped in the 
electromagnetic field 0, El IE & IE or held in a glass 
cell [IS M, [Hill- 
Ultimately, quantum memory should allow to store a 
quantum state of one or several modes of the electromag- 
netic field, and not only states of single photons. Partic- 
ularly promising for this purpose appears to be the off- 
resonant quantum non-demolition (QND) coupling be- 
tween a collective (pseudo)-spin of the atomic ensemble 
and the polarization of light beam [3, EE HE El ■ This 
setting involves an auxiliary coherent laser beam that 
mediates the coupling between an atomic ensemble and 
a light mode whose quantum state should be stored to the 
atoms. Using large atomic ensembles and strong auxil- 
iary laser beams leads to a collective enhancement of the 
atoms-light coupling. Consequently, sufficiently strong 
interaction can be achieved even without a cavity just 
by a single passage of light through the atomic sample 
[lOl [l2l ItH [lSj . which greatly simplifies the experiment. 
The QND coupling has been explored in several landmark 
experiments to entangle states of two atomic ensembles 



[To| , to store quantum state of a light beam onto the col- 
lective spin of atoms [Tp , and to teleport quantum state 
of light onto atoms [19J. Besides quantum information 
processing applications, the QND coupling can be also 
used to generate spin squeezed state of the atomic en- 
semble that can enhance precision of atomic clocks jg, |9J . 

Collective spin degree of freedom of atomic ensemble 
can be described by collective atomic spin operators J x , 
J y and J z . By making the expectation value of one of 
these operators sufficiently large, the other two manifest 
similar algebraic properties as quadrature operators of 
light. The QND coupling between atoms and light can 
then be described by linear input-output transformations 
of the effective quadrature operators of atoms and light 
[IE EE El • This greatly simplifies theoretical analysis of 
atomic memory operation HE HE HE HE H, HI but also 
somewhat restricts possible manipulations with memory. 
Experimentally easily accessible coherent and squeezed 
light beams are described by Gaussian Wigner functions 
and all such states are referred to as Gaussian states. 
Using Gaussian light states and QND coupling, possibly 
combined with homodyne detection on output light and 
feedback, we can implement only Gaussian operations on 
the atomic memory and prepare only Gaussian states in 
the memory starting from initial Gaussian state. 

Certain applications, notably continuous- variable en- 
tanglement distillation, however require non-Gaussian 
operations [H, HE HE HE HI EH Ell. :t is there- 
fore highly desirable to investigate schemes for imple- 
mentation of non-Gaussian operations and filters on the 
quantum-memory state and devise procedures for prepa- 
ration of arbitrary highly nonclassical states of atomic 
quantum memory. Previously, scheme for generation 
of a superposition of two coherent spin states of an 
atomic ensemble has been proposed [33| and a protocol 
for probabilistic noise-free upload of single-photon and 
Schrodinger cat-like states into atomic memory has been 
suggested [33 |. 

In this paper, we propose a scheme for condi- 
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FIG. 1: Atoms- light interaction setup. Light is prepared in 
squeezed vacuum state in the squeezer and then a tiny portion 
of the light beam is reflected from an unbalanced beam split- 
ter BS and impinges on single-photon detector APD. Click of 
APD heralds subtraction of single-photon from the squeezed 
beam. Light in such state accompanied by orthogonally po- 
larized strong coherent beam interacts with atomic ensemble. 
Afterwords, measurement of pl quadrature is performed on 
output light beam using homodyne detector (HD). 



tional preparation of arbitrary coherent superpositions 
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Let Na denotes the number of atoms in the ensemble. 
The quantum memory typically exploits coherence be- 
tween two atomic Zeeman or hyperfine levels | f) and 
||). The Dicke state \n) is then defined as a fully sym- 
metric state of all Na atoms with n atoms in state | f) 
and Na — n atoms in state | J,). In the limit of large Na, 
the states \n) become formally equivalent to the n-photon 
Fock states of light. 

The general idea behind our protocol is to manipu- 
late the atomic state by QND interaction with light that 
has been previously prepared in a specific highly non- 
classical quantum state. This technique allows us to em- 
ploy QND interaction to implement operations on atomic 
ensemble that are not easy to perform directly on atoms, 
but are more feasible on light. In particular, we show 
that by using light beam prepared in photon-subtracted 
squeezed vacuum state [3ol l3ol l37j we could implement 
an operation on the atomic memory that is similar to 
single-photon subtraction/addition. 

This elementary non-Gaussian operation can be com- 
bined with coherent displacements of atomic state and 
repeated several times to conditionally generate a wide 
class of superpositions (fT]). The resulting scheme is anal- 
ogous to protocols for generation of arbitrary superposi- 
tions of Fock states of traveling light beams by repeated 
photon addition or subtraction [3a . [3{| , but it exhibits 
important differences due to the QND coupling and the 
need for conditioning on the outcomes of homodyne de- 
tection on the output light beam. We analyze in detail 
influence of various relevant experimental parameters on 
the performance of the protocol and show how it can be 
optimized in order to maximize the fidelity of generated 
state for a given probability of success. 

The rest of the paper is organized as follows. In Sec- 
tion II we will review the basics of QND coupling between 



light and atoms and derive a general formula for the op- 
eration performed on atoms. In Sec. Ill we explain how 
the elementary non-Gaussian operation can be combined 
with magnetic-field induced atomic displacement opera- 
tor in order to generate various superpositions of states 
|0) and |1). In Section IV, the procedure is generalized 
to multiple repeated application of non-Gaussian opera- 
tion and coherent displacement, which allows to gener- 
ate wide class of superpositions ^ . Results of numerical 
simulations and strategies for optimization of the proto- 
col are presented in Sec. V. In Section VI we show that 
by controlling which atomic quadrature couples to light 
we can in principle generate arbitrary complex superposi- 
tions (UJ) . In Sec. VII we discuss an alternative procedure 
for direct single-step preparation of atomic state by using 
a specifically prepared highly non-classical state of light, 
which is then imprinted onto atoms. Finally, Section VIII 
contains a brief summary and conclusions. 



II. ATOMS-LIGHT INTERACTION 

Our scheme for preparation of highly non-classical 
states of the atomic ensemble is based on the off-resonant 
quantum non-demolition interaction between light beam 
and collective atomic spin. The light beam propagating 
along z axis consists of a strong coherent vertically po- 
larized mode and a horizontally polarized mode prepared 
in a pure quantum state \4>l)- The atoms are prepared 
by optical pumping in a coherent spin state with all spins 
pointing along the x axis such that the x-component of 
the collective atomic spin J attains a macroscopic value, 
(J x ) = FNa, where Na is the total number of atoms in 
the cloud. Under this condition we can replace opera- 
tor J x with its mean value in the commutation relations 

J y , J z ] = iJ x and define the effective atomic quadratures 

1( 



J II 



J, 
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PA 



2(J X ) 



2(Jx) 



satisfying canonical commutation relations 
[xa,Pa] = o- 



(2) 



(3) 



Similarly, the quadratures of the horizontally polarized 
light mode are defined such that [xl,Pl] = i/2 holds. 
With this normalization the wave-function of vacuum 
reads 4> v&c {x) — (2/tt) 1 / 4 exp(— x 2 ), which facilitates fur- 
ther calculations. Using quadrature operators one can 
write down the effective QND interaction Hamiltonian 
between light and atoms in the form of 0, [H, E3] 



-ffQND = hKXLXA, 



(4) 



where R is the interaction constant. 

Suppose that before the interaction atoms and light 
are in the initial pure states \4>a) an( i \4>l), respectively. 
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Light beam passes through the atoms and is subsequently 
subjected to homodyne detection of the Pl quadrature. 
The operation performed on the atomic state correspond- 
ing to a particular measurement outcome pi, can be ex- 
pressed as 



9 = (p L \ exp (-LffQNDVM I0l 



(5) 



where \pl) is the eigenstate of Pl with eigenvalue pl- 
Working in the ^-representation and using the fact that 
(pl\xl) — -l= e ~ l2xLPL we obtain 

e = ^/ c~ 2 ™^e- 2 ™ XAX ^ L {x L )&x Ll (6) 

where k = Rt/2. This formula can be written in a concise 
form 



6 = $(ra A +Pl), 



(7) 



where $(pl) denotes the Fourier transform of the wave 
function 4>l{xl) = {xl\4>l)- 

Formula ((7| reveals that in order to apply a non- 
Gaussian operation on the atomic state we need non- 
Gaussian state of light beam \4>l)- Recently, highly non- 
Gaussian states of traveling light beam exhibiting nega- 
tive Wigner function have been gene rated by photon sub- 
traction from squeezed vacuum [35l l36l W\ . As schemat- 
ically shown in Fig. [TJ squeezed light beam is generated 
in an optical parametric amplifier and then impinges on 
a highly unbalanced beam splitter BS with nearly unit 
transmissivity. Click of the single-photon detector APD 
indicates with high probability a subtraction of single 
photon because probability that two or more photons are 
reflected is negligible. The resulting conditional trans- 
formation can be described by the action of annihilation 
operator on the initial squeezed state, 



S L {r)\0i 



where Sl(t) = exp 



a L S L (r)\0i 



(8) 



^t2 



denotes the squeezing 



2 y^L — u zJ 

operator and |0l) is the vacuum state. The success prob- 
ability of the whole procedure depends on the efficiency 
of the detector and also on the reflectivity of the beam 
splitter. One may increase the success probability of the 
operation by using more reflective beam splitter. This 
will however result in higher probability of reflection of 
two photons and therefore in lower fidelity of the pre- 
pared state. In the limit of very low reflectance of BS, the 
photon-subtracted squeezed state is formally equivalent 
to squeezed single-photon Fock state and its normalized 
wave-function reads 



4>l{xl) 



1/4 



^' 2 2x 



Lexp i 



(9) 



Performing Fourier transformation we obtain the result- 
ing operator 

Bs(pl) = N{x A + Pl/k) exp [-e(x A + p L /n) 2 ] , (10) 
where e = K 2 e~ 2r and J\f stands for normalization con- 
stant, 



N = 2k I - 

7T 



1/4 



-3r/2 



(11) 



The operator 0g is the key essence of our quantum state 
engineering scheme. One can easily verify that action 
of Os(0) on vacuum atomic state \0a) leads to the final 
atomic state whose wave function reads 



(x a \Qs(0)\0a) 



1/4 



Afx A exp [— (e + l)x A ] . (12) 



We thus obtain a highly nonclassical squeezed state |1) 
in the atomic memory [34j. We can see that the ©s(0) 
operator acts simultaneously as a combination of creation 
and annihilation operator and a squeezing operator. As 
shown in next section, this property can be exploited to 
generate (squeezed) superpositions of higher Dicke states. 



III. COMBINATION WITH DISPLACEMENT 
OPERATOR 

Displacement operator D(a) = exp(aa^ — a*a) can be 
easily implemented on atomic ensemble by application of 
magnetic field resulting in a tiny rotation of the collec- 
tive atomic spin. We propose to combine displacement 
operator and the Os operator to generate various super- 
positions of states |0) and |1). In this and the following 
section we will assume conditioning on the measurement 
outcome p^ = corresponding to operation 0s(0). In 
practice, we have to use a finite acceptance window for 
the measurement outcomes in order to achieve finite suc- 
cess probability. The effect of the width of the acceptance 
window on fidelity of generated state will be discussed in 
Sec. V. 

Consider a sequence of coherent displacement followed 
by non-Gaussian operation Os(0) and another displace- 
ment. The entire operation on the atomic ensemble takes 
the form 



D(a 2 )e s (0)D(a 1 ). 



(13) 



The displacement operator D{oi2) can be propagated to 
the right using D{pij)x = (x — aj)D(ctj), where ctj — 
a,- + ibj and au, bj are real. We obtain 
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D(a 2 )Q S (0)D(ax) = e^ a ^\x A - a 2 ) exp [-e{x A - a 2 f] D(at + a 2 ). (14) 

Assuming the atoms are initially in the effective vacuum state |0) the transformation (|13p prepares the atoms in a 
state 

ipA,out(xA) oc {x A - a 2 ) exp [-(e + l)x\ + 2 (a a + (e + l)a 2 ) x A + 2ix(bi + b 2 )] , (15) 



where an unimportant constant has been neglected. 
Wave function ([To*]) represents a squeezed superposition 
of vacuum and |1) state provided that its Gaussian part 
is centered on origin. This can be accomplished simply 
by putting 



h=b 2 = 0, 
ai = — (e + l)a 2 . 



(16) 



By changing the free parameter a 2 , one can easily modify 
relative amplitudes of states |0) and |1) in the superposi- 
tion (|15p as long as they remain real, i.e. without relative 
phase shift. We can thus prepare in the atomic mem- 
ory arbitrary squeezed state S(s)(c 10) + Ci 1 1)) , where 
co,ci are real and s = |log(e + 1). In current experi- 
ments, k < 1 holds. If the light beam is squeezed in the 
amplitude quadrature, r > 0, then e <C 1 and the re- 
sulting squeezing of the state in atomic memory is small 
and could be neglected. However, if the light beam is 
squeezed in phase quadrature, r < 0, then the resulting 
atomic squeezing may become significant. We will see in 
Sec. V that choosing r < may be advantageous since it 
can provide better trade-off between success probability 
of the protocol and fidelity of the generated state. 

The fact that the target state co|0) + ci|l) is squeezed 
is not necessarily an obstacle and, in fact, it can be de- 
sirable in certain cases. For instance, the squeezed state 
|1) very well approximates superposition of two coherent 
states | a) — | — a) , which is often referred to as Schrodinger 
cat-like state and is a valuable resource for quantum in- 
formation processing. The squeezing can be removed 
by the anti-squcezing operation S(—s) performed on the 
atomic state. This can be accomplished by means of 
an auxiliary light beam that interacts with the atomic 
state several times either sequentially [2(| or simultane- 
ously [25]. If e < 1 then one could also pre-squeeze the 
atomic ensemble before state preparation to compensate 
for the effects of squeezing. Initializing the atomic mem- 
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FIG. 2: Schematized multiple interaction scheme for genera- 
tion of superpositions of first N + 1 Dicke states. The proce- 
dure consists of sequence of N non-Gaussian operations com- 
bined with N + 1 coherent displacements of the atomic state, 
as described by Eq. (fl7)) . 



ory in squeezed state ipA,in(xA) oc exp[— (1 — e)x A ] would 
yield the desired superposition of Dicke states without 
any parasitic squeezing. Note that in order to have the 
Gaussian part of the wave-function centered on vacuum, 
we should now set a\ = —a 2 /(l — e). The initial squeez- 
ing of the state of atomic ensemble can be accomplished 
by QND measurement of the pa quadrature followed by 
coherent displacement of the atomic state proportional 
to measurement outcome [10, \vA . 



IV. MULTIPLE INTERACTION SCHEME 

By repetition of the basic sequence fTS]) one is able to 
prepare superpositions of higher Dicke states. Such gen- 
eralized operation acting on initial atomic vacuum state 
reads, 

D(a N+1 )e s {0)D( aN )e s {0) . . . D{a 2 )e s {0)D( ai )\0) 

(17) 

and is schematized in Fig. [2] Similar calculation as that 
leading to Eq. (fl"2"]) gives us the resulting atomic wave 
function. Defining 



JV+l 

E 

j=k 



a i 



(18) 



and decomposing these cumulative displacements into 
real and imaginary parts, au = fife + ibk, the wave func- 
tion reads 



JV+l 



jV+1 



4>a(xa) oc (xa - fife) exp -eNx 2 A + 2eJ2 



a k XA 



fe=2 



x exp 



fc=2 



N+l 

' > Ofe 

fe=2 



^ fife + 2ib\XA — (xa — fii 



% \2 



Inspired by Eq. (|16[) . one can always set 



JV+l 

fii + e fife 

fc=2 



'3 = 0, Vi, 
0. 



(19) 



(20) 



The Gaussian part of the wave function (jTU]) then be- 
comes centered on origin, 

N+l 

iPa(xa) oc Yl ( X A - fife) exp [- (eN + 1) x\] . (21) 
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FIG. 3: (color online) Average state fidelity as a function 
of success probability is plotted for the target state = 
(|0) + |l})/\/2 and parameters k = 0.5 and r = 1. Red full 
line: basic acceptance window, green dotted line: advanced 
acceptance window, blue dashed line: advanced acceptance 
window and feed-back. 



FIG. 4: (color online) Average fidelity is plotted as a function 
of success probability of preparation for target state = 
(|0) + |1))/V5, K = 0.5, and r = 1 (blue dashed line, light 
squeezed in xl) and r = — 1 (red full line, light squeezed in 
Pi)- 



Wave function of any Dicke state |n) is composed of prod- 
uct of a Hermite polynomial of nth order H n (y2x) and 
a Gaussian. Therefore the wave function of finite super- 
position of Dicke states ([T|) can be expressed as a prod- 
uct of a Gaussian and a polynomial with degree equal 
to the highest Dicke state in the target superposition. 
The polynomial part of the wave function ^target can be 
decomposed into a factorized form, 



N 

-x 2 \ " 



^target (x) OC e ~ ^ C - /=- f 

„=o v2"n! 



H n (V2x) 2 yj 

n oc e Wlx — Rj) 



By comparison with the wave function that results 
from multiple interaction scheme ([2Tj) . we see that we 
are able to engineer a squeezed version of the desired 
state (|22p just by setting correctly displacement param- 
eters ay. We need to find roots Rj of the polynomial 
J2n=o °nH n (V%x) /V2™n! and set bj = and 



(ik 



Rn 



a i 





r 


Rk-i - 


- Rk 




f 1 


1 









,N, 



N 



-R\ — e / ,Rk 
fe=i / 



(23) 



Note that this formula must yield real a,j . Our protocol 
is therefore capable of preparation of any superposition 
of Dicke states as long as all roots Rj remain real. In par- 
ticular, we can prepare arbitrary (squeezed) Dicke state 
|n), because all roots of Hermite polynomial H n (x) are 
real. 

Using the parameters given by Eq. (f2"3"]l . the state 
obtained by the preparation procedure differs from the 



desired target state just by squeezing in the xa quadra- 



ture and can be expressed as S(s)\tp, 



target/ ) 



where now 



s = h \og(eN +1). As discussed in the previous section, 
one may employ an anti-squeezing operation that results 
in rescaling of quadrature operator, 



x A 



X A 



VeN+ 1 



(24) 



After such operation the atomic state will become the 
desired superposition of Dicke states |T]). Provided that 
Ne < 1 one could also pre-squeeze the atomic ensemble 
before the state engineering procedure and start the pro- 
tocol from state %1)a{%a) oc [—(1 — Ne)x A \. The protocol 
would then directly yield the required superposition of 
Dicke states provided that 



ffljv+i — Rn 

a k = R k -i - Rk, k = 2 

N 

-Ri 



....,N, 



a i 



1 - Ne 



^Rk- 



(25) 



fc=i 



V. NUMERICAL SIMULATIONS 

The proposal as it has been presented in the previous 
sections assumed conditioning on pr, = 0. In the ex- 
perimental realization one can not post-select only cases 
where pl — but has to use a finite acceptance win- 
dow to get finite non-zero probability of success. The 
simplest option is to post-select cases, when the result 
satisfies \pi\ < rj, where rj is the acceptance threshold. 
With larger rj the overall probability of success increases, 
because the chance of pl falling into the acceptance win- 
dow is higher. However as rj differs from the prepared 
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FIG. 5: (color online) Joint probability of measurement out- 
comes pli and pl2 in the preparation of Dicke state |2), 
K = 0.5, r = — 1 (colormap). The contour lines indicate the 
fidelity of the prepared state for given measurement outcomes. 
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FIG. 6: (color online) Average fidelity as a function of success 
probability in the process of preparation of |2) state using 
K = 0.5, r = — 1. Red full line: basic acceptance window, 
green dotted line: advanced acceptance window, blue dashed 
line: advanced acceptance window and feed-back. 



state p becomes mixed and its fidelity defined as over- 
lap with pure target state decreases. Thus we obtain 
a trade-off between fidelity and probability of successful 
preparation. The acceptance condition \pi\ < 77 is sim- 
ple, but not optimal because generally fidelity is not a 
monotonic function of \pi\ . One can use a better strategy 
in selecting which values of p L are considered successful 
measurement outcomes. Our advanced acceptance win- 
dows are constructed by choosing some fidelity threshold 
and all values of p L for which this threshold is exceeded 
are considered successful outcomes. By changing the fi- 
delity threshold we get the trade-off between fidelity and 
success probability, because the lower the fidelity thresh- 
old is, the bigger the success probability gets, but the 
average fidelity gets lower. To improve the resulting fi- 
delity even more, one may employ a feed-back strategy. 
For each measurement outcome pl it is possible to find 
corresponding displacement imposed on the atomic state 
that maximizes the state fidelity. 

We present results of numerical simulation in Fig. [3J 
where we plot the trade-off between fidelity and suc- 
cess probability in the case of preparation of = 
(|0) + |l))/\/2. In this simulation we have set k = 0.5 
and r = 1. We can see that for this particular target 
state and parameters, the fidelity of the generated state 
drops quickly if the basic acceptance condition \p L \ < 77 is 
used. The advanced acceptance window provides a bet- 
ter trade-off, as clearly visible in Fig. [3] Further signifi- 
cant improvement can be achieved if we make use of the 
feed-back and coherently displace the state depending on 
the measurement outcome. Combination of advanced ac- 
ceptance window and feed-back yields fidelities exceeding 
90% even for success probability of the order of 50%. 

In some situations, it may be advantageous to use light 
state squeezed in the phase quadrature p L , hence r < 0. 
Fig. 2] illustrates this effect by comparing the trade-offs 



between fidelity and success probability obtained for r = 
1 and r = —1, respectively, with all other parameters 
identical. In both cases, we use advanced acceptance 
windows and optimized feed-back. 

The concept of advanced acceptance window can be 
straightforwardly generalized to preparations of super- 
positions of first N + 1 Dicke states, in which case we 
apply the non-Gaussian operation ©s iV-times and ob- 
tain N measurement outcomes Plj- We evaluate the 
fidelity of the generated state for each set of outcomes 
{Plj} and the preparation is considered successful if the 
fidelity exceeds certain threshold. 

We demonstrate this procedure on the example of 
preparation of (squeezed) Dicke state \2). The protocol 
involves two operations ©s interspersed with three dis- 
placements and we have two measurement outcomes p L ,i 
and pl, 2- The probability of detection of a pair of out- 
comes pl.i and p^ j2 is plotted in Fig. [5] as a function of 
Pl,i and pl,2- The same figure also contains contour plot 
of the fidelity of the prepared state as a function of pl.i 
and pl.2- Advanced acceptance windows are represented 
by areas in the plane defined by these fidelity contours. 
Note that when calculating the fidelity we again assumed 
feed-back in form of coherent displacement of the final 
atomic state and optimized the value of the displacement 
for each pair of outcomes. The trade-off between fidelity 
and success probability for the preparation of the Dicke 
state 1 2) is presented in Fig.[6]which contains comparison 
of three preparation strategies: basic acceptance window 
\PL.j\ < t], advanced acceptance window, and advanced 
acceptance window combined with feed-back. 
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VI. COMPLEX COEFFICIENTS 

In previous sections we have described a procedure 
that is capable of preparation of any superposition of 



atomic Dicke states with real Rj, c.f. Eq. (|22[) . In this 
section we propose a generalized scheme that allows us 
to prepare any complex superposition. The necessary 
degrees of freedom are obtained by coupling the light 
beam to different atomic quadrature at each application 
of non-Gaussian operation O5. Let us consider the ro- 
tated atomic quadrature 



xa,, 



xa cos 4> + Pa sin ( 



(26) 



In the experiment, the atoms-light interaction Hamilto- 
nian (H|) can be changed to Hkxa,$Xl by application of 
magnetic field which rotates the atomic spin about the x 
axis. When the light couples to the the rotated atomic 
quadrature, the interaction operator Qs,r/> is given by 



e 



s,0 oc kx a ,4> exp {-ex A rj> ) 



(27) 



where we assume conditioning on pi = 0. Similarly as 
before, we combine the Qs.<j> operator with displacement 
operators D(a.j). Let us now consider only two repeated 
applications of the Qs, 4, operator enveloped by three dis- 
placement operators. The overall operation reads 



£>(a 3 )e S)02 D(a 2 )e s ,^ J D(ai). 



For the sake of simplicity we will assume the limit of 
strong squeezing of the light beam and/or weak atoms- 
light coupling, e — > 0. Under the condition a\+a 2 + a^ = 
the action of operator (f28|) on atomic vacuum state 
produces a complex superposition of three lowest atomic 
Dicke states, 



|Vu)oc|2) + Cl |l) + co|0), 



(29) 



where c\ and cq represent relative amplitudes between 
states |1) and |2), and |0) and |2), respectively. The two 
complex numbers Co and c\ uniquely specify any super- 
position of states |0), |1) and |2) with non-zero contribu- 
tion of the highest Dicke state |2). After some algebra we 
arrive at the following set of equations for the complex 
parameters Cj, 



c = -^={xie- i ^x 2 e- i ^ + e" 2 *), 
v2 

ci = -- Uzie-** 1 + x 2 e~ 1 ^), 
V2 



(30) 



where x 1 = [{a 2 + a^e" 101 ] and x 2 = 5R[a 3 e"^ 2 ]. 
This set of equations can be effectively reduced to equa- 
tion for a single complex parameter z = x 2 exp (— i(f> 2 ), 



(31) 



We can formally solve Eq. (|3ip as quadratic equation in 
z to obtain 

(32) 

Taking the absolute values of both sides of Eq. (f3"2")l we 
derive an equation for \z\, 




2\/2co [ 1- r-rrr 




(33) 

On the left-hand side we have a monotonically growing 
function satisfying lim z ^ oc \z\ = 00. On the right-hand 
side we have a function that approaches asymptotically a 
finite real constant for large values of \z\ while for \z\ — > 1 
it grows to infinity. In the interval (1; 00) the function on 
right-hand side of Eq. (f33|) is continuous. Therefore there 
has to exist an intersection point between the graphs of 
functions on right-hand and left-hand sides. This point 
is thus a positive real solution of the equation (|33|) . 

We have proved that there is always a solution of the 
equation (f3"3"|) . Moreover, this solution can be determined 
analytically because this equation can be transformed to 
a polynomial equation of the fourth degree, 



diV A + d 3 y 3 + d 2 y 2 + d lV + d = 0, 



(34) 



(28) where y = 



and 



d 3 

d 2 
di 

d 



1. 



-4-2| Cl | 2 , 

6 + 4| Cl | 2 - 4|c | 2 + 2V2(c 2 cS + cf c ), 



= -4-2|ci| 



4|c | 2 | Cl | 2 



■8|co| 



-2\f2(c\c*o +cf c ) 
l + 4|c | 4 -4|c | 2 . 



(35) 



Finding the non-negative real root of this polynomial 
gives the exact solution of the whole problem. We have 
confirmed numerically for a large number of randomly 
generated Co and ci that there is always such root. The 
presented method allows us to prepare arbitrary complex 
superpositions of the first three atomic Dicke states and 
overcomes the limitations of the scheme discussed in pre- 
vious Sections. This technique can be straightforwardly 
extended to generation of arbitrary superpositions of the 
first N Dicke states and the effect of finite e can also 
be included. However, the resulting nonlinear equations 
for the coherent displacements au and the phase shifts 
<j)j become highly complicated and can be solved only 
numerically. 



VII. DIRECT MAPPING APPROACH 

In this section we present an alternative strategy for 
preparation of superpositions of atomic Dicke states. 
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Success probability 

FIG. 7: (color online) Trade-off between fidelity and success 
probability in the case of preparation of |2) Dicke state. The 
Fourier approach is indicated by full red line, the original 
sequential approach by blue dashed line, k, = 0.5. 



This approach is based on the preparation of specific 
state of light beam, which is then imprinted into the 
atomic memory. By using appropriate superposition 
of the first N Fock states of light and conditioning on 
Pl = 0, we can prepare the atoms in a single step in 
arbitrary desired (squeezed) superposition of the first N 
Dicke states. 

The wave-function of target atomic state squeezed by 
a factor of y/1 + re 2 can be expressed as 



4>A (xa) 



2(1 + k 2 ) 



1/4 



JV 



y/2 n nl 



-H, 



n (V2(1 + k 2 )xa) 



(36) 



Since this wave-function should be obtained by applying 
the operation (&(kxa) onto the vacuum atomic state, we 
can determine the required wave function of light <fr(p) 
in the p-representation, 



Explicitly, we have 

1 



1/4 



e XA ip A (x A ). 



(37) 



JV 



Hp) oc ' 



V2"n! 



-.H„ 



^2(1 + « 2 ) 



P e 



This wave function describes a superposition of the first 
N + 1 Fock states, 



\</>l) = ^ u m \m). 



(38) 



m=0 



The coefficients u m can be determined by calculating the 
scalar product (m|^i) in the ^-representation. We arrive 



at the following formula for the non- normalized u m , 

N 



E 

xH. 



\/2 m + n m\n\ 



^2±^pl e-*dp. (39) 



As an example, we have determined the light state \4>l) 
that is required for the preparation of the (squeezed) 
atomic Dicke state |2). Explicit calculation yields 



\4> L ) oc |0>- V2(1 + k 2 )|2). 



(40) 



The advantage of the direct single-step preparation is 
that only one interaction between atoms and light is 
needed. Moreover, we can generate an arbitrary atomic 
state by employing an appropriate state of the light 
beam (|38p. The latter can be prepared from initial vac- 
uum or squeezed state by combination of coherent dis- 
placements and single-photon additions or subtractions 
[H, [39[. However, the probability of success of such 
conditional preparation schemes typically decreases ex- 
ponentially with the maximum number of photons N in 
the superposition. In contrast, if we use the sequential 
preparation scheme described in previous sections, then 
we can wait for successful preparation of the photon- 
subtracted squeezed state before we switch on the cou- 
pling between atoms and light. We still have to condition 
on the outcomes of the homodyne detections, so even for 
this scheme the success rate will decrease with growing 
N and fixed fidelity of the generated state. Which strat- 
egy is optimal thus depends on the particular target state 
and other experimental parameters. 

Figure [7] shows the trade-off between fidelity and suc- 
cess probability in the case of preparation of the |2) Dicke 
state. This trade-off is compared with the trade-off for 
the original sequential strategy. Note that the success 
probability represents only the probability of acceptance 
of the measurement outcomes pl of homodyne detector 
and does not include the probability of generation of the 
required states of light. As discussed above, the total cost 
of preparation of the non-Gaussian states of light would 
be lower in case of the sequential strategy than in case 
of the single-step strategy. In the numerical simulations, 
we assumed advanced acceptance window and optimized 
feed-back on the atomic state. Interestingly, the sequen- 
tial strategy outperforms the direct strategy in this case 
as can be seen in Fig. [7] 



VIII. CONCLUSIONS 

In this paper, we have proposed a protocol capable of 
preparation of a wide class of highly non-classical states 
of atomic quantum memory. We have shown that a non- 
Gaussian operation on the state stored in memory can 
be performed using non-Gaussian state of light, QND 
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coupling between atoms and light and conditioning on 
the homodyne detection performed on the output light 
beam. Based on this elementary non-Gaussian opera- 
tion, we have devised a simple strategy for preparation of 
real superpositions of atomic Dicke states by sequential 
application of coherent displacement and non-Gaussian 
operation <d$- Several numerical simulations have been 
performed to verify the functionality of our scheme. We 
have found that the performance of the protocol can be 
significantly improved by judicious choice of the accep- 
tance window for results of homodyne detection p^ and 
by applying feedback displacement on atoms depending 
on pl. We have also shown that using coupling of light to 
arbitrary rotated atomic quadratures allows us to prepare 
any complex superpositions. Finally, we have presented 
a general single-step scheme for preparation of atomic 
Dicke state based on imprinting a state of light beam onto 
quantum memory. We hope that the obtained results 
will stimulate attempts to experimentally generate highly 
non-classical states of atomic memory and that they will 



be found useful in development of advanced quantum in- 
formation processing schemes with atomic memories re- 
quiring non-Gaussian operations. 
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